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Abstract. We consider some classes of Levy processes for which the estimate of 
Krylov and Safonov (as in [BL02]) fails and thus it is not possible to use the stan- 
dard iteration technique to obtain a-priori Holder continuity estimates of harmonic 
functions. Despite the faliure of this method, we obtain some a-priori regularity 
estimates of harmonic functions for these processes. Moreover, we extend results 
from [SSV06] and obtain asymptotic behavior of the Green function and the Levy 
density for a large class of subordinate Brownian motions, where the Laplace 
exponent of the corresponding subordinator is a slowly varying function. 



I. Introduction 

Recently there has been much interest in investigation of the continuity properties 
of harmonic functions with the respect to various non-local operators. An example 
of such operator £ is of the form 

(Cf)(x) = [ (f(x + h)- f(x) - (V/(z), h)l w<1 ) n(x, h) dh (1.1) 

JR d \{0} 

for / G C(M d ) bounded. Here n: M. d x (M. d \ {0}) — > [0, oo) is a measurable function 
satisfying 

Cl \h\- d - a < n{x,h) < c 2 \h\~ d - a , 
for some constants c±, c% > and a G [0, 2]. 



It is known that for a G (0, 2) Holder regularity estimates hold for ^-harmonic 
functions (see [BL02] for a probabilistic and [Sil06] for an analytic approach) . 
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For a = 0, techniques developed so far are not applicable. One of our aims is to 
investigate this case using a probabilistic approach. 

In many cases the operator of the form (1.1) can be understood as the infinitesimal 
generator of a Markov jump process. The kernel n(x, h) can be thought of as the 
measure of intensity of jumps of the process. 

Let us describe the stochastic process we are considering. Let S = (S t : t > 0) be 
a subordinator such that its Laplace exponent defined by 0(A) = — log (Ee _A51 ) 
satisfies 

lim — — — = for any x > (1.2) 

A^+oo 0'(A) 

for some a G [0, 2] . Let B = (B t , F x ) be an independent Brownian motion in M. d and 
define a new process X = (X t , P z ) in M d by X t = B(S t ) . It is called the subordinate 
Brownian motion. 



Example 1: Let S be a subordinator with the Laplace exponent satisfying 

r m i 

aJJ^o A«/ 2 £(A) 

with a G (0,2) and t: (0, oo) — \ (0, oo) that varies slowly at infinity (i.e. for any 
x > 0, — > 1 as A — > +oo). 

Then (1.2) holds (we just use Theorem 1.7.2 in [BGT87] together with the fact that 
0' is decreasing and 0(A) = f Q </>'(t) dt) and 

n(x, h) = n{h) x \h\- d ~ a £{\h\- 2 ), \h\ 0+ , 

which means that ^-d-^ h ^ stays between two positive constants as \h\ — > 0+ . 

Choosing i = 1 we see that the rotationally invariant a-stable process (whose in- 
finitesimal generator is the fractional Laplacian C = — (— A) 2") is included in this 
class. Other choices of i allow us to consider processes which are not invariant under 
time-space scaling. 



The processes described in Example 1 with I = 1 belong to the class of Levy stable 
or, more generally, stable-like Markov jump processes. The potential theory of 
such processes is well investigated (see [BL02, CK03, SV04, BS05, BK05, RSV06, 
KS07, CK08, MimlO, SztlO, KM11]). For example, it is known that harmonic 
functions of such processes satisfy Holder regularity estimates and the scale invariant 
Harnack inequality holds. Also, two sided heat kernel estimates are obtained for 
these processes. 
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Not much is known about harmonic functions in the case when the corresponding 
subordinators belong to 'boundary' cases, i.e. a G {0,2} (see [SSV06, Mimll, 
Miml2]). For the class of geometric stable processes only the non-scale invariant 
Harnack inequality was proved and on-diagonal heat kernel upper estimate is not 
finite (see [SSV06]). 

The following two examples belong to these 'boundary' cases and are covered by our 
approach. 

Example 2: Let 0(A) = log(l + A). The corresponding process X is known as the 
variance gamma process. In (1.2) we have a = . Moreover, it will be proved (see 
Theorem 4.1) that 

n(x,h) = n(h) x \h\- d , \h\ ->• 0+ . (1.3) 

This example can be generalized in various ways. For example, we can take k G N 
and consider 0& = o . . . o 0. Then (see Theorem 4.1 or Subsection 7.1) 

k times 

n k (x,h) =n k (h) x \h\~ d I log ■ ■ ■ log ^ ■ . . . ■ log log ^ ■ log ^ J , \h\ ->• + . 

\ fe— 1 times / 

Another generalization is to consider 0(A) = log(l + A^ 2 ) for some (5 G (0,2]. The 
process X is known as the /3-geometric stable process and the behavior of n is given 
also by (1.3). 

Example 3: Let 0(A) = log ^^ ■ Then in (1.2) we have a = 2 and (see Theorem 
4.1) 

n(x, h) = n(h) x \h\- d - 2 (log ' , \h\ -> + . 

This behavior shows that small jumps of this process have higher intensity than 
small jumps of any stable process. 

A measurable bounded function /: M d — > R is said to be harmonic in an open set 
D C M d if for any relatively compact open set B C B C D 

f(x) = E x f(X TB ) for any x G B, 

where t b = inf{t > 0: X t B} . 

The main theorem is the following regularity result, which covers cases a G [0, 1). 
The novelty of this result rests on the case a = 0. By B r (xo) we denote the open 
ball with center xq G M. d and radius r > . 
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Theorem 1.1. Let S be a subordinator such that its Levy and potential measures 
have decreasing densities. Assume that the Laplace exponent of S satisfies (1.2) with 
a G [0, 1). Let X be the corresponding subordinate Brownian motion and let d>3. 

There is a constant c > such that for any r G (0, |) and any bounded function 
f : R d — > R which is harmonic in B^O), 



\f(x) - f(y)\ < cll/Hoo j' for all x,y G Br(0) . 



Applying Theorem 1.1 to Example 1, we obtain expected Holder regularity esti- 
mates. Within this example the result is new when the scaling is lost, e. g. 

0(A) = At [log(l + A)] 1_ * . 



The situation is more interesting in Example 2, e.g. for the geometric stable process. 
For this process we obtain logarithmic regularity estimates: 

l/W - /(»)! < c||/I|oolog(r- x ) 



logflx - y\ 



It is still unknown whether Holder regularity estimates hold for harmonic functions 
of this process (or, generally, of the processes belonging to the case a = 0). 



Let us explain why known analytic and probabilistic techniques do not work in the 
case a = 0. The main idea in the proof of the a priori Holder estimates of harmonic 
functions relies on the estimate of Krylov and Safonov. 



In probabilistic setting this estimate can be formulated as follows. There is a con- 
stant c > such that for every closed subset A C B r (0) and x G Bl(0) 

V*(TA<T B(M )>Cj£^ (1.4) 

where Ta = ta? is the first hitting time of A and \A\ denotes the Lebesgue measure 
of the set A . 



Performing a computation similar to the one in the proof of Proposition 3.4 in [BL02] 
(see also Lemma 3.4 in [SV04]) we deduce 

> c 



P*(T A < r Br( o)) > 



r 2 0'(r 2 ) 






<p(r 2 ) 


B r (0)\ 
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If a G (0,2), it can be seen that r - x 1 as r — >• 0+ . This gives estimate 

of the form (1.4) and thus the standard Moser's iteration procedure for obtaining 
a-priori Holder regularity estimates of harmonic functions can be applied (see the 
proof of Theorem 4.1 in [BL02] for a probabilistic version). 

The situation is quite different for a = 0. To find a counterexample we will use the 
following result. 

Proposition 1.2. Let S be as subordinator such that its Levy and potential measures 
have decreasing densities and whose Laplace exponent satisfies (1-2) with a G [0, 1). 
Let X be the corresponding subordinate Brownian motion and let d> 3. 

There is a constant c > such that for every r G (0, 1) and x G -Br (0) 
K(X TBi(0) G B r (0) \ Br(0)) < c r ~ 2 ^ 2 2) . 



For a = it will follow that lim r ^-2) ~ = (see (2.10)). Therefore in this case 
(1.4) does not hold, since 

r hm Po(T Br (o)\B s (o) < r Br{0) ) < lim P (I TBr(0) G B r (0) \ Bj( )) = . 



Considering process X in the setting of metric measure spaces (as in [CK08] or 
[BGK09]) the news feature appears. Theorem 4.1 shows that the jumping kernel of 
the process X is of the form n(x, h) = j(\h\) with 

r- 2 0' (r- 2 ) E r Br(0) 



In the case a = the term - J^ r -2^ - becomes significant. This has not yet been 
treated within this framework. 

The latter discussion shows that the question of the continuity of harmonic functions 
becomes interesting even in the case when the kernel n(x,h) is space homogeneous, 
or in other words, in the case of a Levy process. There is no known technique that 
covers this situation in the case of a more general jump process. 

Our technique relies on asymptotic properties of the underlying subordinator. The 
potential density can be analyzed using the de Haan theory of slow variation (see 
[BGT87]). 
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On the other hand, there is no known Tauberian theorem that can be applied to 
obtain asymptotic behavior of the Levy density fi of the subordinator. For this 
purpose we perform asymptotic inversion of the Laplace transform (see Proposition 
3.2) to get 

n(t) xrV(r 2 ), t^o+. 



These techniques allow us to extend results from [SSV06] to much wider class of 
subordinators whose Laplace exponents are logarithmic or, more generally, slowly 
varying functions. 



Although we do not obtain regularity estimates of harmonic functions for cases when 
a G [1,2], it is possible to say something about the behavior of the jumping kernel 
and the Green function. In this sense, the case a = 2 is also new. For example 
the Green function of the process corresponding to the Example 3 above has the 
following behavior: 

G(x,y) x \x - y\ 2 ~ d log(\x - yl' 1 ), \x - y\ ^0 + . 

We may say that such process X is 'between' any stable process and Brownian 
motion. 



Let us briefly comment the technique we are using to prove the regularity result. In 
Section 2 it will be seen that any bounded function / which is harmonic in -B 2r (0) 
can be represented as 



f(x) = / K Br{0) (x,z)f(z) dz, x G B r (0) , 

J WW 

where K Br ^(x, z) is the Poisson kernel of the ball B r (0). 

The following estimate of differences of Poisson kernel is the key to the proof of 
Theorem 1.1: 



\K Br (o)(xuz) - K Br (o)(x 2 ,z)\ < 
for x±,X2 G Bl(0) (see Proposition 5.3). 



-rf 0((M-r)- 2 ) 
4>(\x- 

<f>(\x-y\- 2 ) 



c \z r 1 i 2,' r<z<2r 
\z\ > 2r 



Similar type estimate has been obtained in [SztlO] for stable Levy processes using 
scaling argument and the explicit behavior of the transition density. In our setting 
there are many cases where the behavior of the transition density is not known and 
the scaling argument does not work. Our idea is to establish the following Green 
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function difference estimates: 



\G(xi,y) - G(x 2 ,y)\ < c 



r 




r 



,d 



1 A 



\*-y\ 



r 



) 



for all y G R d and x±,x 2 G" B r (y) (see Proposition 5.1). 

The paper is organized as follows. In Section 2 we introduce all concepts we need 
throughout the paper. Section 3 is devoted to the study of subordinators. We obtain 
asymptotic properties of Levy and potential densities. In Section 4 asymptotical 
properties of the Green function and Levy density of the subordinate Brownian 
motions are obtained. Difference estimates of the Green function and the Poisson 
kernel are the main subject of Section 5. This type of estimates are the main 
ingredient in the proof of the regularity result in Section 6. In Section 7 we apply 
our results to some new examples. 

Notation. For two functions / and g we write / ~ g if f/g converges to 1 and 
/ x g if f/g stays between two positive constants. The n-th derivative of / (if 
exists) is denoted by p n \ 

The logarithm with base e is denoted by log and we introduce the following notation 
for iterated logarithms: log! = log and log fc+1 = log o log fc for fcGN. 

We say that /: R — y R is increasing if s < t implies f(s) < f(t) and analogously for 
a decreasing function. 

The standard Euclidian norm and the standard inner product in R d are denoted 
by | ■ | and (-,-), respectively. By B r (x) = {y G R d : \y — x\ < r} we denote the 
open ball centered at x with radius r > 0. The Gamma function is defined by 



2.1. Levy processes and their potential theory. A stochastic process X = 
(X t :t> 0) with values in R d (d > 1) defined on a probability space (fl, J 7 , P) is said 
to be a Levy process if it has independent and stationary increments, its trajectories 
are P-a.s. right continuous with left limits and f(X Q = 0) = 1 . 



r(p) = J^t^e-'dt for p > 0. 



2. Preliminaries 



The characteristic function of X t is always of the form 

EexpO(£,X t )} = exp {-*$(£)}, 
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where $ is called the characteristic (or Levy) exponent of X. It has the following 
Levy-Khintchine representation 

$(0 = i<7,0 + ^,0+ / {l-e?™ + i(x,t)l w<1} )lL(dx). 

Here 7 G A is a non-negative definite symmetric d x d real matrix and II is a 
measure on R d , called the Levy measure of X, satisfying 

n({0}) = and / (1 A \x\ 2 )U{dx) < 00 . 



The Brownian motion B = (B t : t > 0) in M. d with transition density po(t,x,y) = 
(47rt)~ d / 2 exp | — I is an example of a Levy process with the characteristic ex- 
ponent = |£| 2 . 

A subordinator is a stochastic process S — (S t : t > 0) which is a Levy process in R 
such that St G [0, 00) for every t > 0. In this case it is more convenient to consider 
the Laplace transform of S t : 

Eexp{-AS t } = exp{-t0(A)}, A > 0. 

The function 0: (0, 00) — > (0, 00) is called the Laplace exponent of S and it has the 
following representation 

0(A) = 7 A+ / (l-e~ xt )/i(dt). 

J (0,00) 

Here 7 > and /i is also called the Levy measure of S and it satisfies the following 
integrability condition: // Ooo )(l /\t)/i(dt) < 00. 

The potential measure of the subordinator S is defined by 



U(A) = E 



l{5 t eA} ^ 



for a measurable Ac [0, 00) . 
The Laplace transform of U is then 

CU(X):= [ e- xt U(dt) = -3— . (2.1) 



(0,00) 



0(A) 



Assume that the processes B and S just described are independent. We define a 
new stochastic process X = (X t : t > 0) by X t = B(S t ) and call it the subordinate 
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Brownian motion. It is a Levy process with the characteristic exponent $(£) 
</>(|£| 2 ) and the Levy measure of the form U(dx) = j(\x\) dx with 



j(r) 



(0,00) 



{Aixt)- d ' 2 exp {--) »(dt). 



(2.2) 



The process X has the transition density and it is given by 

12 ' 



p(t,x,y) 



-d/2 



cxp 



[0,00) 



F - y\ 

4s 



F(S t e ds) 



(2.3) 



When X is transient, we can define the Green function of X by 



G(x,y)= p(t,x,y) dt, x,y e R , x ^ y . 
^(0,00) 

The Green function can be considered as the density of the Green measure defined 
by 



l{x t eA} dt 



G(x,A) = E. 
since G(x, A) = J A G(x, y) dy. 

Using (2.3) we can rewrite it as G(x, y) — g(\y — x\) with 



icl measurable , 



(0,00) 



(A7rt)- d / 2 exp {--\ U(dt) 



(2.4) 



Let D C M. d be a bounded open set. We define the process killed upon exiting D 
X D = (X t D : t > 0) by 

Y D _ j X t t <T D 

1 ~\ d t>r D , 
where d is an extra point adjoined to D. 

Using the strong Markov property we can see that the Green measure of X D is 



G D (x,A)=E x 



^{x t D eA} dt 



l{x t eA} dt 



l{Xt6A} dt 



G{x, A) - E X [G(X TD ,A); r D < 00] . 
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Thus in the transient case the Green function of X D can be written as 

G D {x,y) = G(x, y) - E X [G (X TD , y); t d < oo], x,y e D, x^y. 

Since X is, in particular, an isotropic Levy process it follows from [SztOO] that 

F x (X TBr{0) e dB r {0)) = 0, xeB r (0). 

for any r > and x G B r (0) . This allows us to use the Ikeda-Watanabe formula 
(see Theorem 1 in [IW62]): 

F X (X G F) = [ [ G Br{G) (x,y)j(\z-y\)dydz, (2.5) 

J F JB r (0) 

for x G B r (0) and F C 5 r (0)° . 

C 

Defining a function Kbjq) '■ B r (0) x -B r (0) — >■ [0, oo) by 

K Br (o)(x,z)= G Br ( )(x,y)j(\z-y\)dy (2.6) 

the Ikeda-Watanabe formula (2.5) reads 

¥ X (X G F) = [ K Br{0) (x, z) dz . (2.7) 
Jf 

The function i^_B r (o) will be called the Poisson kernel for the ball B r (0) . 

2.2. Bernstein functions and subordinators. A function 0: (0, oo) — > (0, oo) is 

said to be a Bernstein function if G C°°(0, oo) and (— l) n (f)^ < for all n G N. 
Every Bernstein <fi function has the following representation: 

0(A) = 7i + 7 2 A + / (1 - e- xt )fi(dt), (2.8) 

J(0,oo) 

where 71,72 > and fi is a measure on (0, 00) satisfying C 0oo \(l A t)/i(dt) < 00. 

Using the elementary inequality ye -;y < 1 — e _3/ , y > we deduce from (2.8) that 
every Bernstein function satisfies: 

A0'(A) < 0(A) for any A > . (2.9) 

There is a strong connection between subordinators and Bernstein functions. To be 
more precise, is a Bernstein function such that 0(0+) = (i.e. 71 = in (2.8)) if 
and only if it is the Laplace exponent of some subordinator. If 0(0+) > 0, can be 
understood as the Laplace exponent of a subordinator killed with rate 0(0+) (see 
Chapter 3 in [Ber96]). 
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A Bernstein function is a complete Bernstein function if the Levy measure in (2.8) 
has a completely monotone density, i.e. p(dt) = p(t) dt, where p: (0, oo) — > (0, oo), 
p E C°°(0, oo) and (-1) V n) > for any n G N . 



Let us mention some properties of complete Bernstein function (see [SSV10]). The 
composition of two complete Bernstein function is a complete Bernstein function 
and if is a complete Bernstein function, then 0*(A) = is also a complete 
Bernstein function. 



Assume that S is a subordinator with the Laplace exponent <fi and infinite Levy 
measure. Then <ft* is a Bernstein function if and only if the potential measure U 
has a decreasing density u with respect to the Lebesgue measure. Moreover, if v 
denotes the Levy measure of the subordinator with the Laplace exponent <ft*, then 
u(t) = u(t, oo) for any t > . 



2.3. Regular variation. A function /: (0, oo) — > (0, oo) varies regularly (at infin- 
ity) with index p 6 R if 

/(Ax) 

lim = x p for every x > . 

A^+oo /(A) 

If p — 0, then we say that / is slowly varying. Regular (slow) variation at is 
defined similarly. 

If / varies regularly with index p e R, then there exists a slowly varying function I 
so that /(A) = A^(A) . 

Let i be a slowly varying function such that L(X) = J Q A ^ dt exists for all A > . 
Then L is slowly varying, 

lim ^ = +oo (2.10) 



and 



L(Xx)-L(X) , 
hm — — = log x for every x > . 

A^+oo £[X) 



(see Proposition 1.5.9 a and p. 127 in [BGT87]). 
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3. SUBORDINATORS 



Let S = (S t : t > 0) be a subordinator with the Laplace exponent <fi satisfying the 
following conditions: 

(A-l) there is a G [0,2] such that <p' varies regularly at infinity with index f — 1, 
i.e. 

(j)'(Xx) a_i 

hm = x 2 for every x > , 

A^+oo <p'(X) 

(A-2) the Levy measure is infite and has a decreasing density fi , 
( A-3) the potential measure has a decreasing density u . 

When a = 2 we additionaly assume: 

(A-4) A i — y C$@j) var i es regularly at infinity with index — f . 

Remark 3.1. (a) The most important assumption is (A-l) (and (A-4) when a = 
2). Other assumptions hold for a large class of subordinators (e. g. when 
is a complete Bernstein function with infinite Levy measure), 
(b) By Karamata's theorem (see Theorem 1.5.11 in [BGT87]) it follows from 
(A-l) that varies regularly at infinity with index ^ • 

Proposition 3.2. Let a G [0,2) and let S be a subordinator satisfying (A-l) and 
(A-2). Then 

Proof. Let e > 0. By a change of variable 

POO 

0(A + £)-0(A)= / {e- xt - e- (A+e) ')Mt) dt 



o 



= A 

Since \x is decreasing, the following holds 



/•oo 

1 / e-'(l -e- £X ~ H )fi(\-H)dt. (3.1) 
Jo 



0(A + e) -0(A) > A-V(A _1 ) / e-\l-e- eX ~ H )dt. 



o 



Now we can apply Fatou lemma to deduce 



= hm ^ (A + £) - 0(A) > A-V(A~ l ) fte-Ut 
= \- 2 /i(\- 1 )(l-2e- 1 ). 
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By setting A = t^ 1 we get the upper bound 

H(t) < - ^ J for every t>0. (3.2) 

Now we prove the lower bound. Using (3.1), for any r G (0, 1) and e > we can 
write 

<f>(\ + e) - 0(A) = h + h (3.3) 

with 

h = A -1 f e"*(l - e- £X ~ H )^\-H) dt 
Jo 

POO 

I 2 = A" 1 / e -*(l - e- £A_1 *)/i(A- 1 t) . 

r 

Since // is decreasing, the dominated convergence theorem yields 

limsup — < \~ 2 {i{\- l r) / te"*rft 

e^0+ £ J r 

= (r + l)e- r A- 2 /i(A"V) . (3.4) 



To handle I\ first we use the theorem of Potter (see Theorem 1.5.6 (iii) in [BGT87]) 
to conclude that there are constants c\ > and 5 > such that 

V \ < Cit s for all A > 1 and t<\. (3.5) 
0'(A) 

Therefore, by (3.2), (3.5) and the dominated convergence theorem 

limsup — < limsup / e~* — : t~ 1 6 / (\t~ 1 )dt 

<0'(A) r -^ ZT £t 5 - 1 e-^t. (3.6) 

Combining (3.3), (3.4) and (3.6) we deduce 

0'(A) < <t>'{\) - C \ - x ^ t 5 - l e- 1 dt + (r + l)e" r A- 2 /i(A-V) . 
1 — ^ e Jo 

Furthermore, by choosing r G (0, 1) so that x _^e~ x Jq t^e' 1 dt < ~ we get 

7* 

/i(A-V) > 6 , rr AV(A) for every A > 1 . 
2(r + 1) 
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By (A-l) we can find to £ (0,r) so that ^fer^ > f° r an Y * £ (0, to)- The lower 



bound nwo follows: 



1+- r 

tit) > A( i^N ^'Vft" 1 ) for every t G (0,t ) 
4(r + 1) 



□ 



Remark 3.3. The precise asymptotical behavior of \x when a G (0, 2) can be obtained 
by Karamata's Tauberian theorem. This is not the case when a = 0. 
Under an additional assumption 

t i — y ta at tit) is monotone on (0, T) for some a > and T > 0, 

it is possible to prove the following precise asymptotics in the case a = : 

tit) ~ rV(t _1 ), t ->■ 0+ . 

Proposition 3.4. Let a G [0,2) and let S be a subordinator satisfying (A-l) and 
(A-3). Then 

u(t) ~ — , r- ' \ NO , t ->• o + . 

w r(i-f) ^(t- 1 ) 2 

Remark 3.5. It can be proved that 

i rVfr 1 ) 
u r(i-f) ^(t- 1 ) 2 

similarly as in Proposition 3.2. It is enough to note 

POO 

ip(X + e) - ip(X) = / (e- xt - e- {x+£)t )u(t) dt 
Jo 

With V(A) = ~^y. 

The main reason why we need precise asymptotics of u is to be able to handle the 
case a = 2 by duality. 



Proof of Proposition 3.4- Let us first consider the case a = 0. In this case £(X) = 
X(p'(X) varies slowly (at infinity) and thus it follows from Subsection 2.3 that 0(A) = 
io ^ ^ a l so varies slowly and 

<f>(Xx) - <f>(\) , 
lim — — = log x for every x > . 

A^+oo X(f)'(X) 
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This and (2.1) imply 

cu(^)-cu{{) H i)- H ^4, { i x) 
iHi)H\f " Mi) *(i)" logI ' A ^° + 

for any x > 0. Now we can apply de Haan's Tauberian theorem (see -version of 
[BGT87, Theorem 3.9.1]) to deduce 

U(Xx) - U(X) 

— — H -> logx, A ->■ + . 

If we apply de Haan's monotone density theorem (see [BGT87, Theorem 3.6.8]) we 
finally obtain 

, x ry (r 1 ) 

u(t) - t v ' t 0+ . 

The case a G (0, 2) is already known. We give the proof for the sake of completeness 
and adapt the result to the formula obtained in the case a = . 

Since 

CU(X) 1 



0(A) 

varies regularly at infinity with index — 2-, Karamata's Tauberian theorem (see The- 
orem 1.7.1 in [BGT87]) implies 

1 1 

U([0,t]) ? \ ~r~? — TT ; * -> + • 

u J; r(i-f)0(t- 1 )' 

Then by Karamata's monotone density theorem (see Theorem 1.7.2 in [BGT87]) we 
deduce 



a 1 1 t- 2 (h'(t- 1 ) 

2r(i-f)t0(t- 1 ) r(i-f) ^(t- 1 ) 2 

□ 



Now we consider the case a = 2. 

Proposition 3.6. Lei a = 2 and let S be a subordinator satisfying (A-l), (A-2) 
and (A~4)- Then 

n(t) ~ r 2 (^(r 1 ) - <l>\t~ x )) , t -»■ o + . 
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Proof. Since potential density exists by (A-3) we see that it follows that </>*(A) = 
defines the Laplace exponent of a (possibly killed) subordinator, which we de- 
note by T (see Section 2). 

Note that the subordinator T corresponds to the case of a = 0. If we denote 
potential density of T by v, then 

v(t) = fi(t, oo), t > 0. 

Proposition 3.4 yields 

Ks) ds ~ ^i} )2 \ t + . (3.7) 

By assumption (A-4) we know that t h-> (0*)'(t _1 ) varies regularly at with index 
1 and thus t >->■ * ' varies regularly at with index —1. 

Now we change variable in the integral on the left-hand side in (3.7) and conclude 



This gives (r = t x ): 

cis r 2 (4>*)'(r 
s 2 0*(r) 2 



/ /i(s _1 ) — ~ — VY \ M y , r — t- oo . 
Jo 



.r — )• oo 



Note that the right-hand side is now regularly varying at infinity with index 1 and 
thus by Karamata's monotone density theorem (see Theorem 1.7.2 in [BGT87]) we 
deduce 

^(r -1 ) r(0*)'(r) 
r 2 0*(r) 2 

Going back (t = r _1 ) we conclude 

, > r* (</>*)' (r 1 ) 

'^° + 



□ 



Proposition 3.7. Let a = 2 and Zet S be a subordinator satisfying (A-l) and (A-3). 
Then the following is true 
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Proof. By (2.1) we get 

cu{x) = W) ~ WW' x -* +0 ° 

and thus by Karamata's Tauberian theorem (see Theorem 1.7.1 in [BGT87]) it 
follows that 

since A i— > A0'(A) varies regualrly at infinity with index 1. By applying Karamata's 
monotone density (see Theorem 1.7.2 in [BGT87]) theorem we deduce 

«(*)~ 5? P).*->o+. 

□ 



4. Levy density and Green function 



Let S be a subordinator as in Section 3 and let X be the corresponding subordinate 
Brownian motion in M. d with d > 3. Our aim is to establish asymptotical behavior 
of the Levy density and Green function of X . 

Recall that the Levy density of X is of the form where j is given by (2.2). 

Theorem 4.1. Assume that S satisfies (A-l) with some a G [0,2] and (A-2). If 
a e [0,2), then 

j(r) x r~ d -y (r~ 2 ), r ->■ + . 
If a = 2 and (A-4) holds, then 

j(r) x r-- d - 2 {r 2 (p{r-- 2 ) - 0'(r" 2 )) , r ->■ + . 

Proof. This result follows directly from Proposition 3.2 and Proposition 3.6 together 
with Lemma A.l, where o = j, b — 1 + f . 

For a G [0,2) the slowly varying function is given by i(t) = t^~ l <p'{t~ l ). When 
a = 2 we take 



£{t)=tcj>{t )-0(t ) = -—__, 
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which varies slowly, since </>*(A) = is slowly varying by (A-l) and Karamata's 
Theorem (see Theorem 1.5.11 in [BGT87]) and has a derivative that varies regularly 
with index 1 by (A-4). □ 

The Green function of X is of the form G(x,y) = g(\y — x\), where g is given by 
(2.4). 

Theorem 4.2. Assume that S satisfies (A-l) with some a G [0,2] and (A-3). Let 
d>3. Ifae [0,2), then 

9(r) x r + . 



If a = 2, then 



9(r)xr " +2 ^ xr "^- r ^ + 



Proof. We use Lemma A.l with a — j, b — 1 — ~, Proposition 3.4 and Proposition 
3.7. 



When a G [0, 2) we define £(t) = 1 - which varies slowly at . 

In the case a = 2, we let £(t) = ^pry or £{t) = t ^-i^ which both vary slowly at 
0. " □ 

Using asymptotical results from this section we can now prove the proposition that 
gives a counterexample for the estimate of the Krylov and Safonov. 



Proof of Proposition 1.2. By (2.5), 

F x (X TBr(o) eB r (0)\Br(0))= [ f G B (x,y)j(\z-y\)dy 

* ' JB r (0)\Br(0) JBr(O) ^ 

T T 

= h + h- 
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Using Theorems 4.1 and 4.2 it follows that 
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h= / G Br (x,y)j(\z-y\)dydz 

JB r (0)\Br(0) JB 3r (0) 3 
3 TS 



<j(^)\B r (0)\Br(0)\ / g(\y\)dydz 

JB r (0) 

~ 1 d>(r~ 2 ) 



On the other hand, 



G B (x,y)j(\z - y\) dydz 



B r (0)\Br(0) JBr{0)\B 3r (Q) 



<s(h)/ / j(\y\)dydz. (4.1) 

JB r (0)\B^ (0) JBr( Z ) 

To estimate the inner integral, note that B^(z) C -Bi(O) \ B\ z \_r(0) for any z G 



B r (0) \ Sr(0) and so, by Theorem 4.1, 



Br(z) J\Z\- 



j(\y\) dy<c 2 s- 3 ^'(s- 2 ) ds < c 3 cf>((\z\ - r j)~ 2 ) . (4.2) 



Thus, by Theorem 4.2 and (4.1) 



4 



r~ 2 <f)'(r~ 2 ) 
" C5 0(r- 2 ) • 



In the last equality we have used Karamata's theorem (see Theorem 1.5.11 in 
[BGT87]) and the fact that a G [0, 1) . □ 
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5. Difference estimates 

Let X be the stochastic process in M. d as in Section 4 and assume that d > 3. In 
particular X is transient. 

In this section we prove the difference estimates of the Green function and the 
Poisson kernel. 

Although we are slightly abusing notation, we set G(x) := G(0, x) = g(\x\). 
Proposition 5.1. There is a constant c > such that for every r e (0, 1) 

\G(x) - G{y)\ < cg{r) (l A for all x,y £ B r (0) . 



Proof. Assume first that \x — y\ < |. By the mean value theorem it follows that for 
any t > there exists $ = y, t) G [0, 1] such that 



it — p it 



< 



2/ 



F - 2/1 



where in the last line the following elementary inequality was used 

se < 2e 2 , s > . 
Then \x + — x)| > \x\ — d\y — x\ > ~ implies 



_M1 _\yV 
e 4t — e it 



By (5.r 



\G{x) - G{y)\ < {^y d ' 2 / r d / 2 

./o 

< 2(47r)- d/2 |a;-|/| 



e « — e 4t 



u{t) dt 



(5.1) 



Since u is non- increasing and varies regularly at with index ^ — 1, by Lemma A.l 
we see that there is a constant c\ > so that 

r d/2 ~ 1/2 e~^u(t) dt < cir- d+1 u(r 2 ) for every r e (0, 1) . 



\G(x)-G(y)\<c 2 g(r)^. 



Theorem 4.2 yields 
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When \x — y\ > |, 

\G(x)-G(y)\<G(x) + G(y)<2g(r) 

since \x\, \y\ > r . □ 

Proposition 5.2. There is a constant c > sitc/j that for all R £ (0, 1), r £ (0, -|], 
yeB R (0)andx 1 ,x 2 eB 4 (0)\B r (y) 



|G Bji (o) (xi, ?/) - G Bjl (o) (x 2 , y) | < c#(r) ( 1 A 



Proof. By symmetry of the Green function, 

G BR (o)(xi,y) = G Ba{0) (y,Xi) = G(x t - y) - Ey[G(X TBaW - x { )] 
= G(x i -y)-E y [G(X TBR(0) -x i )}, 

for i £ {1,2}. Now the result follows from Proposition 5.1. □ 

Proposition 5.3. There is a constant c > such that for any r £ (0, 1) and 

x,yeBr(0) : 

(i) if z £ B 2r (0) \ B r (0), then 

^((N-r)- 2 ). 



\K Br(0) (x,z) - K Br[0) (y, z)\ <c\ 
(ii) if z £" -B 2r (0), i/ien 



Proof. In the estimate 

, v) I j(\z — v\) dv 



\K Br{0 )(x,z) - K Br{0) (y,z)\ < / \G Br (p)(x,v) - G Br{0) (y, 

J B r (0) 

we split the integral into three parts: 



h= j \G Br (o)(x,v) - G Br (o)(y,v)\j(\z - v\) dv 

•> B 2 \ X -y\{x) 



h = I \G Br (o)(x,v) - G Br (o)(y,v)\ j(\z- v\) dv 

>Br(x)\B 2 \ x _ yl {x) 



B r (0)\Br(x) 



\G Br ( )(x,v) - G Br ( )(y,v)\j(\z -v\)dv. 
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For the first part we obtain 



h < / G Br (o)(x,v)j(\z-v\)dv+ G Br{0) (y,v)j(\z-v\)dv 

J B 2 \ x - y \{x) J B s \ x _ y \(y) 

<2j(#)/ G{v)dv< Cl ^^- v (5.2) 

for any z B r (0). We have used Theorem 4.2 to get the last inequality in (5.2). 
In order to estimate I2 we split the integral in the following way. We let N = 



log 2 



and write 



N f 

h < ^2 / \G Br (o)(x,v) -G Br (o)(y,v)\ j(\z-v\)dv. 

n=l ^ B 2" + i\ x -y\( X )\ B 2^\ X - V \( X ) 

Now, for each n G {1, . . . , N} we can apply Proposition 5.2 (with the corresponding 
radii (2 n — l)\x — y\ and r) to get 

/ \G Br ( 0) (x,v) - G Br ( 0) (y,v)\j(\z-v\)dv 

JB nn+1 , ,(x)\Bon\„_„\(x) 



V+l |as— w | ( x )\ B 2 n \x-y\ ( x ) 

< Co 

2™ - 1 



g((2 n -l)\x-y\) , 

^ c 3 ™ — 1 / j(k-^l)^. 

2" + 1 \x-y\( X ) 



By Theorem 4.2 



g ((2--l)|,- i/ |) <c4 ,((2--l)|,- y |) „ e{li2i ... iW}i 

»(f-2/I) ^(k — 2/1) 

with r?(r) = ^~ d ~ 2 ^- 2 )2 • 

Noting that 77 varies regularly at zero with index a — d < 0, the uniform convergence 
theorem for regularly varying functions (see Theorem 1.5.2 in [BGT87]) gives 

^((2 W - l)\x -y\) < c ^ 2n _ for aU n e N and , _ - < 1 
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By Theorem 4.2 and (2.9) g(\x — y\) < ,,, ^ ,_ 2) and so 



<t>{\x-y\- 2 ) 

N 



h < c 6 - lr-'-'gdx - y\)(2 n+1 \x - y\) d j ( 



n=l 

N 



< r , -A 2 / V O^- 1 ) 71 
" 7 0(I^-1/|- 2 )^ 



It remains to estimate 13. Applying Theorem 5.2 we get 
h<c 8 g(r) 1 -?—^ [ j{\v\)dv 

r JB r (z) 

\x - y\<p (\x - y\- 2 ) r~ d f 

^M^nL m)dv - (5 ' 3) 

In the last inequality we have used the theorem of Potter (cf. [BGT87, Theorem 
1.5.6 (hi)]) to conclude that for 5 < 1 — a there is a constant As > such that 

\x - y\<f>(\x - yn ^ A ( \x-y\\ M < 



r0(r~ 2 ) 

since r 1— >■ r0(r~ 2 ) varies regularly at zero with index 1 — a . 
Since 

j(\v\)>j (^) for all v E B r (z) and z e B 2r (0) c 
it follows from (5.3) that 



<p(\x-y\ 2 ) 

On the other hand, for z G B 2r (0) \ B r (0) we deduce from 

B r {z) C ,63(0) \ B {zhr (0) 

(similarly as in (4.2)) that 

j(\v\)dv<c 12 <f>((\z\-r)- 2 ) 

B r (z) 
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By (5.3) 



h < ei 3 \z\- d ^P f or all z G B 2r (0) \ B r (0) . 

<f>{\x-y\ 2 ) 



□ 



6. Regularity of harmonic functions 



Recall that (2.7) gives the representation for any bounded function /: ~R d — > R that 
is harmonic in i? 2r (^o) : 



f(x) = E x flX. 



B r (x ) 



K Br (x )(x, z)f{z) dz, x e B r (x Q ) • (6.1) 



Proof of Theorem 1.1. By (6.1) 

\f(x)-f(y)\< 



K 



B 2r (0) 



Bar (0)1 



x,z) -K B2r (o){y, z)\ dz. (6.2) 



It remains to estimate the integral in (6.2), which we split in the following way 



h 
I 
I 



B 4r (0)\B 2r (0) 



\ K B 2r (o)(x,z) - K B2r(0) (y,z)\ dz 



2= / \K B2r{0) (x,z) - K B2r{0) (y,z)\ dz 

Bl(0)\S 4 r(0) 



= / \ K B 2r (o) (x, z) - K B2r{0) (y, z) | dz 

Bi(0) c 



In order to estimate I x we use Proposition 5.3 (i). More precisely, 



h < 



Cl 



\x - y\ 
c 2 



\z\- d <j) ({\z\ -2rY 2 ) dz 



< 



< 



\x-y\ 

c-2 



B 4r (0)\B 2r (0) 

rV((*-2r)" 2 ) dt 



2r 



\x - y\ 

c 3 



2r 



—Mr)- 1 <j>(r 2 )dt 



\x - y\ 



where in the last inequality we have used Karamata's theorem (see the 0-version of 
Theorem 1.5.11 in [BGT87]). 
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We estimate J 2 and J 3 with the help of Proposition 5.3 (ii). Since the Levy measure 
is finite away from the origin, 

/ 3 < n Q , a / i($)dz< . ( . C5 , 2y 

<P{\x-y\ 2 ) J Bl (or v 1 <t>{\x-y\ 2 ) 

Also, 

where in the last inequality we have used Theorem 4.1. □ 



7. Examples 

In this section is to illustrate our results by some examples. 



7.1. (Iterated) Geometric stable processes. This class of examples belongs to 
the case of a = 0. 

Let (3 G (0,2]. We define a family of functions {</>„: (0, oo) — > (0, oo) : n G N} 
recursively by 

0i(A) =log(l + A /3/2 ), A > 
n+ i = 0i o n , n G N. 

The function 0! is a complete Bernsetin function. Since complete Bernstein functions 
are closed under operation of composition, (p n belongs to this class for every n G N. 

Let 5* n be a subordinator with the Laplace exponent cj) n . S 1 is known as the geomet- 
ric —stable subordinator. We call S n the iterated geometric f-stable subordinator. 
The corresponding subordinate Brownian motions X n will be called (iterated) geo- 
metric /3-stable processes. 

As already remarked in [SSV06], these processes show quite different behavior com- 
pared to the one of stable processes. Our contribution to this class of examples is 
that now we can obtain behavior of the Levy density as a special case of Theorem 
4.1 (even for iterated geometric stable processes). 
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The Levy density of X n is comparable to 

1 1 



\x 



MogfcOl l ) 



which is almost integrable. We can say that (intially) this process jumps slower than 
any stable processes. 

This can be also seen from the behavior of the Green function: 

1 n_1 1 

Gtx,y)x- 5— — TT 7 — - as \x — y\ — > + . 

\x-y\ d log n (\x- y I" 1 ) *± log k (\x - y\ x ) 

As a consequence, ^o T B r (o) x log h=r\ as r — » 0+. Therefore X n needs (on average) 
more time to exit ball B r (0) than any stable process or Brownian motion. 

Theorem 1.1 implies the following a-priori local regularity estimates of harmonic 
functions: 

\f(x) - f(y)\ < cll/lloolog^r- 1 ) , 1 for all x,y € Br(0) 

l og n {\x-y\ L ) 

and any bounded function / which is harmonic in -B r (0) . 

This tells us that the modulus of continuity is bounded by a logarithmic term. It is 
still an open problem whether these harmonic functions satify a-priori local Holder 
continuity estimates. 



7.2. Conjugates of (iterated) geometric stable processes. This class of ex- 
amples corresponds to the case a = 2. 

Let ip n (\) = ^ 7^ , where n are as in Subsection 7.1. 

Since 4> n are complete Bernstein functions, ip n are also complete Bernstein func- 
tions. Therefore, there exist (killed) subordinators T n with the Laplace exponent 
ip n . Killing will not affect the behavior of the Levy and potential densities of T n 
near zero. 

In this case the Levy density of the corresponding subordinate Brownian motion Y n 
behaves near the origin as 

1 1 



|xr 2 log^M- 1 ) fMog^M- 1 
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Note that the integrability conditions of the Levy measure are barely satisfied in 
this case. 

Comparing this behavior to the behavior of the small jumps of the a-stable process, 
we see that small jumps of Y n are more intensive. 

Another interesting feature of this process is the following behavior of the Green 
function: 

G(x,y)^\x-y\ 2 ~ d \og n (\x-y\- 1 ) as \x - y\ ->■ + . 

In this sense the process Y n is 'between' stable processes and Brownian motion, 
since their Green functions are given by 

= Ca \x-y\ a - d and G {2) = c a \x - y\ 2 ~ d . 



Appendix A. Asymptotical properties 

In the appendix we prove a technical lemma which is used throughout the paper. 
Lemma A.l. Let w: (0, oo) — > (0, oo) be a decreasing function satisfying 

w{t) x rt(t), t ->■ 0+, 
for a function £: (0, oo) — > (0, oo) that varies slowly at and b > . 

If p > 1 and a > 0, then 

POD 

I(r)= r p e~^w(t)dt, r > 0, 
Jo 

satisfies 

I(r) x a' p - b+1 r- p+1 w{r), r -> + . 

Proof. Change variables yields 

I{r) = (ar)~ p+1 J™ e -H p - 2 w (y) dt (A.l) 
By assumptions, there are constants ci, c 2 > and r > such that 

IV ( (XT' ) 

c x ar h < -. / < c 2 gT & for every r G (0,r ) . (A.2) 
w(r) 
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Let us first prove the upper bound. Using the fact that w is decreasing, by (A.l) 
and (A. 2) we get 

I(r) < {ar)- p+1 J e- l t p - 2 w (ar) dt + (ar)- p+1 e'H^w (y) dt 
< c 2 (ar)- p+1 a~ b w(r) J e~H p - 2 dt + J e -H p+b ~ 2 dt 

a -p-b+l r -p+l 

for every r G (0, r ) . 



< c' 2 a- p -°^r- p+1 w(r) 



The lower bound follows similarly: 

/oo /*oo 
e-H p - 2 w (ar) dt > Cl (ar)- p+l a- b w(r) / e' 1 ^ 2 dt 

= c[a- p - b+1 r- p+1 w(r) 
for every r G (0, r ) . 



□ 
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